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ABSTRACT 



GRB050904 is very interesting since it is by far the most distant GRB event 
known to date {z = 6.29). It was reported that during the prompt high energy 
CN ■ emission phase, a very bright optical flare was detected, and it was temporal co- 

^ ' incident with an X-ray flare. Here we use two models to explain the optical flare, 

•n . One is the "late internal shock model" , in which the optical flare is produced by 

the synchrotron radiation of the electrons accelerated by the late internal shock, 
^ and the X-ray flare is produced by the synchrotron-self-Compton mechanism. 

O ■ The other is the external forward-reverse shock model, in which the optical flare 

is from the reverse shock emission and the X-ray flare is attributed to the central 
engine activity. We show that with proper parameters, a bright optical flare can 
appear in both models. We think the "late internal shock model" is more favored 
c/3 ■ since in this model the optical flash and the X-ray flare have the same origin, 

which provides a natural explanation of the temporal coincidence of them. In the 



forward-reverse shock scenario, fits to the optical flare and the late afterglow sug- 
gests that the physical parameters of the reverse shock are much different from 
that of forward shock, as found in modeling the optical flash of GRB 990123 
previously. 

Subject headings: Gamma Rays: bursts— ISM: jets and outflows-radiation mech- 
anisms: nonthermal 



1. Introduction 

Gamma-Ray Bursts (GRBs) are bright flashes of high energy photons usually lasting 
about several seconds. They are by far the most luminous objects in the universe, they 
emit so large amounts of energy (up to 10^^ ergs) and thus can be detected up to very high 
redshifts (z > 5). 
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GRB050904 was detected by the Burst Alert Telescope (BAT) onboard Swift on 2005 
September 4 at 01:51:44 UT (Cummings ct al. 2005). It was a long (< 500 seconds duration 
in BAT), multi-peaked, bight burst, the 15 - 150 keV fiuence was (5.4 ±0.2) x 10^^ erg cm~^, 
the spectrum can be described by a power law with a photon index ~ —1.34, its redshift 
has been measured by several groups (Haislip et al. 2005; Antonelli et al. 2005; Price et al. 
2005), z = 6.29 makes it to be by far the most distant GRB discovered to date. 

Boer et al. (2005) reported that they detected a very bright optical flare during the 
prompt high energy emission phase, and at the same time there is an X-ray flare. It is 
widely believed that the reverse shock synchrotron radiation usually peaks in the optical/lR 
band, and this emission has been successful in interpreting the early optical emission from 
GRB990123 (Akerlof et al. 1999; Sari & Piran 1999; Wang et al. 2000; Fan et al. 2002; 
Zhang et al. 2003; Nakar & Piran 2005), GRB021211 (Fox et al. 2003; Li et al. 2003; Wei 
2003; Kumar & Panaitescu 2003), GRB041219a (Blake et al. 2005; Fan et al. 2005b) and 
GRB 050525a (Blustin ct al. 2005; Shao & Dai 2005). However in this reverse shock model, 
it is expected that the emission has a negligible contribution in the X-ray band (However, see 
Fan & Wei 2005). Strong optical flare accompanying an X-ray flare may also be accounted 
for by the "late internal shock model" (Fan & Wei 2005). Originally, that model has been 
proposed to interpret the X-ray flare detected in GRB 011121 (Piro et al. 2005) and many 
XRT X-ray flares (Burrows et al. 2005; Zhang et al. 2005; Nousek et al. 2005). 

The optical afterglow light curve of GRB050904 cannot be described by a simple power 
law, between ~ 3 hours and 0.5 day after the burst, the fading of the afterglow can be 
described by a power law with index -1.36. However after this time the light curve flattened 
to a temporal index of -0.82 (Haishp et al. 2005). Tagliaferri et al. (2005) have found a 
break in the light curve at time tb ^ 2.6 day, which may be the jet effect. In this Letter, 
we try to explain the optical flare with two models, i.e. the reverse shock emission and the 
late internal shock model, and then we will fit the afterglow light curve including energy 
injection and jet effects. 

2. Explanation of the optical flare 

2.1. The late internal shock model 

In the standard shock scenario, the prompt gamma-ray emission is produced by the 
internal shock, and the burst duration is determined by the active timescale of the central 
engine. However, some authors suggest that the activity of the central engine may be much 
longer than the GRB duration, which can give rise to some signatures in multi-wavelength 
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afterglows (Dai & Lu 1998; Zhang & Meszaros 2001; Granot et al. 2003; loka et al. 2005). 
Furthermore, it has been proposed that the Fe hne observed in some GRB X-ray afterglows 
are produced by the late time energy injection (Rees & Meszaros 2000; Gao & Wei 2005). 

Fan & Wei (2005) first proposed the late internal shock model to account for the bright 
X-ray flares detected in many GRBs. Here we will show that the late internal shock model 
not only can produce the X-ray flare, but also can produce the optical flare with proper 
parameters. 

Following Fan & Wei (2005), the typical synchrotron radiation frequency can be esti- 
mated by 



8.5 X io^^(^)^4^(r.. - rfi'^'I^Lll, 

V-^Htl^ Hz, (1) 



where L-m is the outflow luminosity, V sh is the Lorcntz factor of the internal shock, F is 
the Lorcntz factor of the emitting shell, bt is the observed typical variability timcscale, 
and ee arc the energy fractions occupied by the magnetic field and electrons, respectively. 
Here the convention = has been adopted in cgs units throughout the text. 

The coohng Lorentz factor is 7e,c ^ 7.7 x 10^(1 + z)/[{l + Y)TB'^5t], where Y = [-1 + 
^^/Y^^4x€^JeB]/2 is the Compton parameter, x ~ min{l, (^'m/i^c)*-^"^'*^^} (Sari & Esin 2001). 
Then the cooling frequency is 



1.6 X 10^^ 

rl6ti(i 



7.29 ' 

Y)-^ Hz 



'e-/^[r..(r,,-i)]-3/^L;3,^ 



(2) 



The synchrotron-self-absorption frequency is about (Li & Song 2004; Fan & Wei 2005) 



» 2.9 X io"(i±i)-/'4/'iir.»(r.» - 

Caar.-'"«r'" Hz (3) 

where Lgyn is the synchrotron radiation luminosity. The maximum flux of synchrotron 
radiation is F^nax ~ 3-\/3$p(l -|- z)NemeC^crT^B/{327r'^qeDl), where qe is the charge of elec- 
tron, A^e — LjnSt/[{l + z)rmpc'^ is the total number of emitting electrons, $p is a function of 
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p, for p — 2.5, $p ~ 0.6 (Wijers Sz Galama 1999). is the luminosity distance, we adopt 

(Qm,^a,^) = (0.3,0.7,0.71). Then for the case Uc < < ^ohs < the observed flux at 
frequency fobs should be 

^Z,29.3<^^f (l + ^-'mJy (4) 

Now we turn to the observation. Boer et al. (2005) reported that they detected a bright 
optical flare at frequency Vobs = 3 x 10^"^ Hz, the peak flux is 48 mJy. Meanwhile, the Swift 
XRT data shows that there is also a peak in the X-ray light curve at nearly the same time 
with the optical flare, which suggests that the optical flare and the X-ray peak may have the 
same origin. The slope of the X-ray spectrum is about -1/2, and the flux at IKeV is about 
0.08 mJy. 

In our late internal shock model, if we take the values as follows: = 0.4, es = 0.02, 
Lm = 10^2 ergs^S T = 200, T,/, = 1.6, 6t = 20s, then we find z/^ ~ 6.3 x 10^^ Hz, 
Ua ~ 8.4 X 10^^ Hz, t'c ~ 1.2 X lO^^Hz, so it is in the fast cooling phase, between and 
the spectrum takes the form F^, oc and at the observed frequency (3 x 10^^ Hz) 

the flux is 49 mJy, which is quite consistent with the observation. In addition, with the 
values of eg and e^, the Compton parameter y ~ 4, then the synchrotron photons will be 
Compton scattered to high energy, the energy spectrum between 10^^ Hz and 10^^ Hz is 
also F^, oc z/~^/^, and we can estimate the flux at IKeV to be about 0.06 mJy, which is also 
consistent with the observation well. 



2.2. The reverse shock model 

After the internal shock phase, as the fireball is decelerated by the circumburst medium, 
usually a pair of shocks develop (Meszaros & Rees 1997; Sari & Piran 1999; Kobayashi 2000). 
The early optical afterglow lightcurve is usually composed of the contributions from both 
the forward (FS) and the reverse shocks (RS). With that model, the very early optical/IR 
flash following GRB 990123, GRB 021211, GRB 041219a and GRB 050525a could be well 
modelled by assuming the physical parameters are quite different for the FS and RS (Fan 
et al. 2002; Zhang et al. 2003; Kumar & Panaitescu 2003; McMahon et al. 2004; Fan et 
al. 2005b; Blustin et al. 2005). For example. Fan et al. (2002) performed a detailed fit to 
the optical flash of GRB 990123 data and obtained el = 4.7e^ = 0.6 and = 400e^ = 0.4, 
where the superscripts "r" and "f" represent RS and FS, respectively. 
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Boer et al. (2005) found that both the optical flare and the gamma-ray burst of GRB 
050904 were as energetic as those of GRB 990123 (in the rest frame of the GRBs). If other 
parameters (including the initial Lorentz factor of the ejecta, the number density of the 
interstellar medium n) arc similar for these two events, then the resulted shock parameters 
should be similar, too! So it is very likely that in the current case the shock parameters of 
FS and RS are also different. 

Recently, Yan et al. (2005) have developed a code to calculate the GRB afterglow 
light curves, including the FS and the RS emission components. In the current calculation, 
there arc two novel effects have been taken into account. One is that in previous works, the 
Lorentz factor of the outflow as well as the comoving number density of particle are assumed 
to be constant. This may not be the case since in the standard fireball model, the gamma- 
ray burst is from the internal shocks. The detected gamma-ray lightcurve is so variable 
that the involved outflow may be variable, too (both the Lorentz factor and the particle 
number density). In order to model the optical plateau (Boer et al. 2005) and partly for 
convenience, in this work we assume the outflow can be approximated as two parts. Their 
bulk Lorentz factors, isotropic energies and widths are (^/(i), -Z?jso(i), ^(i)) and (?7(2), Eiso{2), 
A(2)), respectively. The other is the more reliable calculation of the arriving time of the RS 
emission. We take the emitting time of the flrst 7— ray photon as our zero point of time. On 
the line of sight (6* = 0), a gamma-ray photon jp arriving us at t implies that the distance 
of corresponding electron (i.e., point P, at which the bulk Lorentz factor is 77) to the initial 
outflow front is ~ ct/{l + z). The radial distance of the FS front to the central engine is 
Rp when the RS crosses point P. At that time, the width between photon 7p and point P is 
~ (1 — (3r^)Rp, where = ^/T^-TJrf■ . Therefore, the arriving time of the RS emission from 
point P should be t + (1 + 2;)(1 — j3rj)Rp/c. It is straightforward to extend this calculation to 
the cases of 7^ 0. It is found that the /—band flare of GRB 050904 can be well reproduced 
with the following parameters (see the insert of Fig. 1): the isotropic energy of the outflow 
is r](i),2 = 380, Eiso{i),54 = 0.4, A(i),i2 = 1.3, r](2),2 = 800, Eiso{2),5A = 0.3, A(2),i2 = 0.7, 
n — 3 cm~^, el = 0.6, and = 0.4. It is surprising to see that the resulting reverse shock 
parameters are nearly the same as those of GRB 990123 (Fan et al. 2002). 

Can the X-ray flare be from the RS, too? The answer is negative. Firstly, as shown 
by Fan & Wei (2005), the dechne of the X-ray emission of the RS can not be steeper than 
^-(2+p/2)^ which is inconsistent with the observation. Secondly, now the reverse shock region 
is significantly magnetized, so the RS emission in X-ray band should also be dominated by 
the synchrotron radiation. Thus the X-ray band emission should be an extension of the 
optical emission. However, the observation shows the optical-to- X-ray bands emission can 
not be described by a simple synchrotron spectra (Boer et al. 2005). Therefore, the X-ray 
flare accompanying the optical flare should be attributed to the activity of the central engine 
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in the RS model. 



3. Fits to the late J— band afterglow 

The multi-wavelength afterglow light curves (especially the J— band one) of GRB050904 
have been detected (Haishp et al. 2005; Taghaferri et al. 2005 and the references therein). 
Between ~ 3 hours and 0.5 day after the burst, the fading of the afterglow can be described 
by a power law with index -1.36. After that time the light curve flattened to a temporal index 
of -0.82. A break appears at time ti ~ 2.6 day, which suggests that the outflow may be a jet. 
In this Letter we pay more attention to the optical flattening. We note that at the observer 
time t ~ 10^ — 10^ s, there arc strong X-ray flares (Cusumano et al. 2005; Price ct al. 2005; 
Watson et al. 2005). Fan & Wei (2005) suggested that when the moderate rclativistic outflow 
powering the X-ray flare catched up with the initial GRB ejecta, a flattening would occur 
in the long- wavelength afterglow light curve. In the calculation, we assume that between 
i ~ 4 X 10^ seconds and f ~ 1.5 day, significant part of energy has been injected into the 
decelerating GRB ejecta. Similar to Zhang et al. (2005), the energy injection rate has been 
taken as dEinj/[dt/{l + z)] — Ac^ {t / to)~^'^ , where A is a constant. We take A — when 
there is no energy injection. With the energy injection, the equation (8) of Huang et al. 
(2000) should be replaced by 

^ _ (1 - l^)dm + A(t/to)-°-5[dt/(l + z)] 
^ Mej + em + 2(1 - e)-fm ' ^ ^ 

where 7 is the bulk Lorentz factor of the GRB ejecta, M^j is the rest mass of the initial GRB 
ejecta, m is the mass of the medium swept by the GRB ejecta (which is governed by dm = 
AnR'^nmpdR, where rup is the rest mass of proton, dR — 7(7 -|- 1/7^ — l)cdt/{l + z), e = xse 
is the radiation efficiency. With the dynamical evolution of the ejecta, it is straightforward 
to calculate its FS emission (e.g., Huang et al. 2000; Yan et al. 2005). 

The fits to the J— band data (taken from Haislip et al. [2005] and Tagliaferri et al. 
[2005]) are presented in Fig. 1. It is found that the data can be well modeled with the 
following parameters: -Biso,54 = 0.7, n = 3 cm~^, = 0.15, = 0.001, A = 7 x 10^^ ergs s"-^, 
^0 = 4 X 10^ s, and the jet angle 9j — 0.054. Note that the value of 9j is obtained from 
fitting the afterglow light curve, not from the simple analytic relation. Comparing with the 
reverse shock parameters derived in §2.2, the shock parameters of the FS and the RS are 
quite different, as that found in GRB 990123 (Fan ct al. 2002; see also Zhang et al. 2003). 
The isotropic energy of the 7— rays is ~ 5 x 10^'^ ergs and the derived 6j = 0.054, so the 
geometry corrected energy should be ~ 7 x 10^° ergs, which is typical for the GRBs detected 
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by BeppoSAX, HETE-2 and Swift. In our treatment, the flattening is caused by the late 
time energy injection. The total isotropic energy injected into the GRB ejecta is ~ 6 x 10^^ 
ergs. 



4. Discussion and conclusion 

The bright optical flare has been detected in GRB050904, which is as bright as the 
optical flash of GRB990123 (in the rest frame of bursts) and seems to be accompanied by 

an X-ray fiarc (Boer et al. 2005). Here we explored two possible models to account for that 
observation. One is the "late internal shock model", in which the optical flare is produced 
by the synchrotron radiation of the electrons accelerated by the late internal shock, and the 
X-ray fiarc is produced by the synchrotron-self-Compton process. ^ The other is the external 
forward-reverse shock model, in which the optical flare is from the reverse shock emission 
and the X-ray flare is attributed to the central engine activity. We show that with proper 
parameters, a bright optical flare can appear in both models. 

In the forward-reverse shock scenario and with late time energy injection, we have 
modeled the optical flare as well as the late J— band afterglow numerically. The resulted 
shock parameters of the forward/reverse shocks are = 4eg = 0.6 and = 400e^ = 0.4, 
respectively. They are quite similar to those found in GRB 990123 (Panaitescu & Kumar 
2001; Fan et al. 2002), which is a natural result in view of the similarity between these two 
GRBs and their optical flares (in the rest frame of bursts). 

As for the reverse shock emission, previous works usually assumed that the physical 
parameters are uniform, this greatly simplify the calculation. But in reality, the observed 
gamma-ray emission light curve is much variable, so it is very likely that the involved outflow 
should be also variable. We notice that if the parameters are uniform, then before the peak 
time the flux rises quickly, cannot account for the observed plateau (Kobayashi 2000; Boer 
et al. 2005). Here just for simplicity, we divide the outflow into two parts. We expect that 
in the realistic case, the outflow should be non-uniform, so the parameters should have a 
continuous distribution within the shell, but the calculation is complicated. 

Although both the "late internal shock model" and the reverse shock emission can 
account for the observed optical flash and the X-ray flare. We think the "late internal shock 



^However, in some cases the synchrotron emission may peak in kcV energy band, then the IC component 
would peak at GeV energy band (unless the outflow is highly magnetized, as suggested by Fan, Zhang & 
Proga (2005a)), which may be detectable for the upcoming GLAST. This possibility will be discussed in 
great detail elsewhere. 
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Fig. 1. — Modeling the /-band flare (the insert) and the J-band afterglow of GRB 050904. 
In the insert, the /—band flare data (the filled rectangles) are taken from Boer et al. (2005), 
the dashed line is the theoretical light curve of the reverse shock emission. The J— band 
afterglow data (the filled circles) are taken from Haislip et al. (2005) and Tagliaferri et al. 
(2005). The solid line is the theoretical light curve of the J— band afterglow. 
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model" is more favored since in this model the optical flash and the X-ray flare have the 
same origin, which provides a natural explanation of the temporal coincidence of them. In 
the late internal shock model, it is needed that after the prompt 7^ray burst phase, the 
central engine could re-start. Recently two models have been proposed for the production 
of late energy injection (King et al. 2005; Pcrna ct al. 2005; MacFadyen et al. 2005 
proposed another model to account for the X-ray flare in short GRBs). While in the reverse 
shock model, the temporal coincidence of the optical flash and the X-ray flare can only be 
regarded as fortuitous. In addition, we note that in the late internal shock model, the typical 
synchrotron radiation frequency strongly depends on the parameters, such as F, Fgh, -^m, 
etc., and for different burst sources it is natural that these parameters are different, so we 
expect that the late internal shock model not only can produce the optical or X-ray flare, 
but also can produce the flare at other wavelength, such as at the ultraviolet or infrared. 
Meanwhile we predict that the synchrotron-self-Compton process may produce emission at 
high energy band (~ GeV). 

Despite its high redshift, the optical afterglow of GRB050904 is not peculiar with respect 
to other GRBs. Recently Zhang et al. (2005) and Nousek et al. (2005) analyzed the X-ray 
afterglows of many GRBs, they found some features (the X-ray flares, the flattening of the 
light curve, a late time break) occurred in a good fraction of GRBs. These features are 
consistent with the afterglow of GRB050904. We suggest that the progenitor of GRB050904 
may be not quite different from that of other GRBs in view of these similarity. 

We thank the referee for her/his helpful comments. This work is supported by the 
National Natural Science Foundation (grants 10225314 and 10233010) of China, and the 
National 973 Project on Fundamental Researches of China (NKBRSF G19990754). 
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